FRACTIONAL CALCULUS OF VARIATIONS 
FOR A COMBINED CAPUTO DERIVATIVE 
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Abstract. We generalize the fractional Caputo derivative to the frac- 
tional derivative '^D"''', which is a convex combination of the left Ca- 
puto fractional derivative of order a and the right Caputo fractional 
derivative of order /3. The fractional variational problems under our 
consideration are formulated in terms of D^'^ . The Euler-Lagrange 
equations for the basic and isoperimetric problems, as well as transver- 
sality conditions, are proved. 
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1. Introduction 

The history of Fractional Calculus (FC) goes back more than three 
centuries, when in 1695 the derivative of order a = 1/2 was described by 
Leibniz. Since then, the new theory turned out to be very attractive to 
mathematicians and many different forms of fractional operators were in- 
troduced: the Grunwald-Letnikov, Riemann-Liouville, Riesz, and Caputo 
fractional derivatives (see, e.g., [MllMlESj), and the more recent notions 
of [25], [l5], and [22l[23]. Besides mathematics, fractional derivatives and 
integrals appear in physics, mechanics, engineering, elasticity, dynamics, 
control theory, economics, biology, chemistry, etc. (cf. [H |T3l [30l [36] and 
references therein). The FC is nowadays covered by several books (e.g., 
[2T| [24l [29l [32| [Ml [39]) and a large number of relevant papers (see, e.g., 

[a El [m [191 [221 [211 [211 [sg [n] ) . 



This is a preprint of a paper whose final and definite form has been published in: Tract. 
Calc. Appl. Anal., Vol. 14, No 4 (2011), pp. 523-537; DOI: 10.2478/sl3540-011-0032-6. 
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The calculus of variations is an old branch of optimization theory that 
has many applications both in physics and geometry (cf. [28l ESj |42l |43] 
and references therein). Apart from a few examples known since ancient 
times such as Queen Dido's problem (reported in The Aeneid by Virgil), 
the problem of finding optimal curves and surfaces has been posed first 
by physicists such as Newton, Huygens, and Galileo. Their contemporary 
mathematicians, starting with the Bernoulli brothers and Leibniz, followed 
by Euler and Lagrange, invented the calculus of variations of a functional 
in order to solve those problems. Fractional Calculus of Variations (FCV) 
unifies the calculus of variations and the fractional calculus, by inserting 
fractional derivatives into the variational integrals. This occurs naturally 
in many problems of physics or mechanics, in order to provide more accu- 
rate models of physical phenomena (see, e.g., [HEllQ]). The FCV started 
in 1996 with the work of [37]. Riewe formulated the problem of the cal- 
culus of variations with fractional derivatives and obtained the respective 
Euler-Lagrange equations, combining both conservative and nonconserva- 
tive cases. Nowadays the FCV is a subject under strong research. Differ- 
ent definitions for fractional derivatives and integrals are used, depending 
on the purpose under study. Investigations cover problems depending on 
Riemann-Liouville fractional derivatives (see, e.g., [H EJ El [l2l \l7\ [T8]). 
the Caputo fractional derivative (see, e.g., [21 UHl EZ] ) , the symmetric frac- 
tional derivative (see, e.g., [25]), the Jumarie fractional derivative (see, e.g., 
[a [3 [221 Eg EI]), and others (see, e.g., [31 HSl [HI [20] ) . Although the liter- 
ature of FCV is already vast, much remains to be done. 

In this paper we extend the notion of the Caputo fractional derivative 
to the fractional derivative '^-D°'^, which is a convex combination of the 
left Caputo fractional derivative of order a and the right Caputo fractional 
derivative of order /3. This idea goes back at least as far as [25], where 
based on the Riemann-Liouville fractional derivatives the symmetric frac- 
tional derivative was introduced. Klimek's approach is obtained choosing 
our parameter 7 to be 1/2. Although the symmetric fractional derivative of 
Riemann-Liouville introduced by Klimek is a useful tool in the description 
of some nonconservative models, this type of differentiation does not seems 
suitable for all kinds of variational problems. The hypothesis that admissi- 
ble trajectories y have continuous symmetric fractional derivatives implies 
that y{a) = y{b) = (cf. [38]). Therefore, the advantage of the fractional 

derivative '^D"'^ via Caputo lies in the fact that using this derivative we 
can describe a more general class of variational problems. It is also worth 
pointing out that our fractional derivative '-"D"'^ allows to generalize the 
results presented in |2]. 
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The text is organized as follows. Section [2] presents some preliminar- 
ies. In Section [3] we introduce the fractional derivative and provide 
the necessary concepts and results needed in the sequel. Our main results 
are stated and proved in Section HI The fractional variational problems 
under our consideration are formulated in terms of the fractional deriva- 
tive ^ D'^'^ . We discuss the fundamental concepts of a variational calculus 
such as the Euler-Lagrange equations for the basic (Subsection 14. ip and 
isoperimetric (Subsection l4.3p problems, as well as transversality conditions 
(Subsection 14. 2p . We end with conclusions in Section [5l 



In this section we review the necessary definitions and facts from frac- 
tional calculus. For more on the subject we refer the reader to [Ml ESI [Ml 
[39j. Let / G Li{[a, b]) and < q < 1. We begin with the left and the right 
Riemann-Liouville Fractional Integrals (RLFI) of order a of a function /. 
The left RLFI is defined by 



where r(-) represents the Gamma function. Moreover, alxf — xlj^f = / if 
/ is a continuous function. 

Let / € AC{[a,b]), where AC{[a,b]) represents the space of absolutely 
continuous functions on [a,b]. Using equations ()2.ip and (|2.2p . we define 
the left and the right Riemann-Liouville and Caputo derivatives as follows. 
The left Riemann-Liouville Fractional Derivative (RLFD) is given by 



2. Preliminaries 




(2.1) 



while the right RLFI is given by 




(2.2) 




(2.3) 



X G [a, b] 



the right RLFD by 




(2.4) 




X G [a, b] 
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the left Caputo Fractional Derivative (CFD) is defined by 



= jI "^fi^)^ ^ e K^]' 

and the right CFD by 

^^D^fix) = -— ^ [\t - xy^^jm 

r(l - a) dt 



(2.5) 



X ■ 



(2.6) 



where a is the order of the derivative. The operators (j2.1|) - (j2.6|) are ob- 
viously linear. We now present the rule of fractional integration by parts 
for RLFI (see for instance [H]). Let < a < 1, p > 1, g > 1, and 
l/p + l/q<l + a. If g e Lp{[a, b]) and / G Lq{[a, b]), then 

rb r-b 

g{x)J^f{x)dx = / f{x)J^g{x)dx. (2.7) 



In the discussion to follow, we will also need the following formulae for 
fractional integrations by parts: 



b . rb 

1-a 



g{x)^,D''J{x)dx= f{x)^ll^-^g{x)\lZ + / f{x):,D'^g{x)dx, 

2.8) 



f g{x)'iDU{x)dx= -f{x)jl-'^g{x)\lJ^+ f f{x)aD^g{x)dx. 

J a J a 

They can be easily derived using equations (|2.3p - (|2.6p . the identity 
(|2.7|) . and performing integration by parts. 

3. The fractional D'^'^ operator 

Let a, P € (0, 1) and 7 G [0, 1]. Motivated by the diamond-alpha opera- 
tor used in time scales (see, e.g., [261110]), we define the fractional derivative 
operator ^ D'^'^ by 

= 7?^: + (1 - iM , (3.9) 
which acts on / S AC {[a, b]) in the following way: 

^<^/(x) = j^D'^Jix) + (1 - 7)^^f /(x). 

Note that 

^D^'^m = ^xD^Jix), ^D^'^fix) = ^D:fix). 
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For f = [/i,...,/iv] : [a,b] ^ M^, (iV G N) and /, G AC{[a,b]), 
i = 1, . . . , N, we have 

The operator (j3.9p is obviously linear. Using equations (|2.8p it is easy 
to derive the fohowing rule of fractional integration by parts for 



+ (1-7) -f{x)Jl~''g{x) 



x=a 

x=b 



x=a 



+ f f{x)D{'^^g{x)dx, (3.10) 

J a 

where -Cfl^ = (1 — j)aDx + j xD^, which acts on / as 

Dtyix) = (1 - 7)al)f /(x) + 7 .^?/(x) • 

Let D denote the set of all functions y : [a, b] such that '^D^'^y 

exists and is continuous on the interval [a,b]. We endow D with the fol- 
lowing norm: 

||y||i,oo := max ||y(x)|| + max \fD"'^y{x)\\, 

a<x<b a<x<b 

where || • || stands for a norm in R^. 

Along the work we denote by diK, i = 1, . . . ,M (M € N), the partial 
derivative of function K : M.^^ — ?> M with respect to its ith. argument. 

Let A G M*". For simplicity of notation we introduce the operators [■] 
and {■}x defined by 



[y](x) = (x,y(x),^Z);'''y(x)) , 
{y}\ix) = (x,yix),^'D'^'^y{x),Xi,, 



4. Calculus of variations via '^D^ 

We are concerned with the problem of finding minima (or maxima) of 
a functional : P — > M, where 2? is a subset of D. The formulation of a 
problem of calculus of variations requires two steps: the specification of a 
performance criterion and the statement of physical constraints that should 
be satisfied. 
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A performance criterion J ^ also called cost functional (or cost), must 
be specified for evaluating the performance of a system quantitatively. We 
consider the following cost: 

J{y)= [\[y]ix)dx, 
J a 

where x € [a,h] is the independent variable, often called time; y(x) € M'^ 
is a real vector variable, the functions y are generally called trajectories or 
curves; '-^ D^'^y{x) G stands for the fractional derivative of y{x); and 
L e C^{[a,, h] X M2^;M) is called a Lagrangian. 

Enforcing constraints in the optimization problem reduces the set of 
candidate functions and leads to the following definition. 

Definition 4.1. A trajectory y G D is said to be an admissible 
trajectory provided that it satisfies all of the constraints along interval 
[a, 6]. The set of admissible trajectories is defined as 

^ {y € D : y is admissible}. 

We now define what is meant by a minimizer of on X>. 

Definition 4.2. A trajectory y G P is said to be a local minimizer 
for JT" on P if there exists 5 > such that J{y) < Jij) for all y G P with 

||y-y||i,(X) < 5. 

The concept of variation of a functional is central to the solution of 
problems of the calculus of variations. 

Definition 4.3. The first variation of J" at y G D in the direction 
h G D is defined as 

e^o e de 

provided the limit exists. 

Definition 4.4. A direction h G D, h 7^ 0, is said to be an admissible 
variation at y G 2? for JT" if 

(i) SJ'{y; h) exists; and 

(ii) y + eh G X> for all sufficiently small e. 
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The following well known result (see, e.g., |42t Proposition 5.5]) offers a 
necessary optimality condition for problems of calculus of variations based 
on the concept of variations. 

Theorem 4.1 ([42j). Let J be a functional defined on V. Suppose 
that y is a ]oca7 minimizer for J on T). Then, 5J{y, h) = for eacii 
admissible variation h in y. 

4.1. Elementary problem of the '-^ D"^'^ calculus of variations. Let 

us begin with the following problem: 

J{y)= I L[y]{x)dx — ^ min (4.11) 

J a 

over all y S D satisfying the boundary conditions 

y{a)=y\ y{b)=y\ y°,y^GM^. (4.12) 

Next theorem gives the fractional Euler-Lagrange equation for the 
problem (l4in) - (i4T2] l. 

Theorem 4.2. Let y = {yi, . . . , i/n) be a local minimizer to problem 
(I4.1ip - ()4.12p . Then, y satisfies the system of N fractional Euler-Lagrange 
equations 

a,L[y](x)+i?ff;ajv+iL[y](x) =0, i = 2,...iV + l, (4.13) 
for all X € [a, b]. 

Proof. Suppose that y is a local minimizer for J'. Let h be an arbi- 
trary admissible variation for problem (|4.1ip - (14.12p . i.e., hi{a) = hi(b) = 0, 
i = 1,...,N. Based on the differentiability properties of L and Theo- 
rem 14. 1^ a necessary condition for y to be a local minimizer is given by 



= 0; 

£ = 



that is, 



'N+l 

_i(x) 

.i=2 

N+1 



+ dN+iL[y]ixfD'^'^hi.i{, 



i=2 



dx = 0. (4.14) 
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Using formula (I3.10p for integration by parts in the second term of the 
integrand function, we get 



■Af+l 



1=2 



hi^i{x)dx 



+ 7 



'N+l 

E 

.i=2 



hi_i{x)xll "5Ar+iL[y](x) 



x=b 



(1-7) 



E 

.i=2 



hi^i{x)all '^dN+iL[y]{x) 



x=b 



0. (4.15) 



Since /ij(a) = hi{b) = 0, i = 1, . . . , N , hy the fundamental lemma of the 
calculus of variations we deduce that 



diL[y]{x) + D^,:^^dN+iL[y]{x) 



0, 



2,...,N + 1, 



for all X € [a,b] 



□ 



Observe that if a and /3 go to 1, then a^x a-^x can be replaced 
with and ^L>f and ^D^ with (see, e.g., [34]). Thus, if 7 = 1 or 
7 = 0, then for a, /3 — )• 1 we obtain a corresponding result in the classical 
context of the calculus of variations (see, e.g., [42]). 

4.2. ^ D^'^ transversality conditions. Let Z € {1, . . . ,N}. Assume that 



y(a) = y'', yi{b) = y° 



1, . . . , N, i I, but yi{b) is free. Then, hi{b) is 



free and by equations (|4.13|) and (|4.15|) we obtain 

7./^"5,+iL[y](x) -{l-j)all-^dN+i+iL[y]i, 



x=b 



Let us consider now the case when y(o) = y", yi{b) = y\ 



(4.16) 
= l,...,iV, 

i I, and yi{b) is free but restricted by a terminal condition yi{b) < y\. 
Then, in the optimal solution y we have two possible types of outcome: 
yi{b) < y\ or yi{b) = y\. If y{b) < y\, then there are admissible neighboring 
paths with terminal value both above and below yi{b), so that hi{b) can 
take either sign. Therefore, the transversality conditions is 



^jl-^di^iL[y]{x) -{l-^)Jl-Pdn+i+im{x) 







x=b 



(4.17) 



for yi{b) < y^. The other outcome yi{b) = y\ only admits the neighbor- 
ing paths with terminal value yi{b) < yi{b). Assuming, without loss of 
generality, that hi{b) > 0, this means that e < 0. Hence, the transversal- 
ity condition, which has it root in the first order condition (|4.14p . must 
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be changed to the inequahty. For a minimization problem, the < type of 
inequahty is cahed for, and we obtain 



7: 



4-^di+iL[y]{x) -(1 - ^)al',-''dN+i+iL[y]{x) 



x=b 



< 



(4.18) 



for yi{b) = yl*. Combining ()4.17p and ()4.18p . we may write the following 
transversality condition for a minimization problem: 



x=b 



yiib) < yl 

imib) - y\) \lJl-''di+,L[y]{x) -(1 - -f)Jl-PdN+i+iL[y]{x) 



<o. 



x=b 



0. 



4.3. The ^ D"^'^ isoperimetric problem. Let us consider now the isoperi- 
metric problem that consists of minimizing (j4.1ip over all y G D satisfying 
r isoperimetric constraints 



g\y)= f &[y]{x)dx = l^, j = l, 

J a 



(4.19) 



where € C^{[a,b] x j = l,...,r, and boundary conditions 

()4.12p . Necessary optimality conditions for isoperimetric problems can be 
obtained by the following general theorem (see, e.g., \A2\ Theorem 5.16]). 

Theorem 4.3. Let J , ■ ■ ■ ,Q'^ be functionals defined in a neigh- 
borhood of y and having continuous first variations in this neighborhood. 
Suppose that y is a local minimizer of (j4.1ip subject to the boundary con- 
ditions ()4.12p and the isoperimetric constrains ()4.19p . Then, either: 
(i) for aU y € D, j = 1, . . . , r, 



6g\yM) sgHy-M) 

6g\y;h^) 6g\y;h^) 
6g''-{yM) 6g' {y;h^) 



5g\y;hn 

6gHT,i^n 

6g'- iy;W) 



(4.20) 



or 



(a) there exist constants Aj G M, j = 1, . . . , r, for which 



5J{y; h)=J2 ^jSg'iy, h) Vh G D. 
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Note that condition (ii) of Theorem 14.31 can be written in the form 

5 (^J{y; h) - \,g\y- h) j = 0. (4.21) 

Suppose now that assumptions of Theorem 14.31 hold but condition (i) 
does not hold. Then, equation (|4.2ip is fulfilled for every h G D. Let 
us consider function h such that h(a) = h(6) = and denote by F the 
functional 

r 

Then, we have 
= 57-(y;h) = ^^T{y + eh%=^ 



'N+l N+l 

diF{y}x{x)hi.i{x) + dN+iF{y}x{xfD''/hi^i{a 



i=2 1=2 

)2N mr 



dx, 



where the function F : [a, b] X M^^^ X M^' ^ R is defined by 



F{yUx) = L[y]{x)-Y^jG^y]{x). 

i=i 

On account of the above, and similarly in spirit to the proof of Theorem l4.2l 
we obtain 

diF{y}x{x)+D^^:^^dN+iF{y}x{x)=0, i = 2,...N + l. (4.22) 

Therefore, we have the following necessary optimality condition for the 
isoperimetric problem: 

Theorem 4.4. Let assumptions of Theorem 14.31 iioid. If y is a lo- 
cal minimizer to the isoperimetric problem (j4.1ip . (j4.12p and (j4.19p . and 
condition (|4.20p does not hold, then y satisGes the system of N fractional 
Euler-Lagrange equations (j4.22p for all x € [a, b]. 

Suppose now that constraints (|4.19p are characterized by inequalities 

gHy)= f G^[y]{x)dx<lj, j = i,...,r. 

J a 

In this case we can set 

{G'[y]{x) - ^) dx + J\cp,{x)fdx = 0, 
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j = 1, . . . , r, where (pj have the same continuity properties as Ui. Therefore, 
we obtain the fohowing problem: 

rb 



Jiy) 



L{x,y{x),^ D"''^y{x),(f){x)) dx 



mm 



where 



[(x), . . . , (j)r{x)], subject to r isoperimetric constraints 



G^[y]ix) 



I. 



dx = 0, j = l,...,r, 



and boundary conditions (j4.12p . Assuming that assumptions of Theo- 
rem 14.41 are satisfied, we conclude that there exist constants Xj £ M, 
j = 1, . . . , r, for which the system of equations 

D^'^^dN+iF{x,y{x),'^ D'^''^y{x),Xi, . . .,\r,4>{x)) 

+ aiF(x,y(x),^Z)"'''y(x),Ai,...,A„0(x)) =0, (4.23) 

i = 2, . . . , iV + 1, where F = L + Y:"j=i XjiG^ + <P]) and 

Aj</),(x) = 0, i = l,...,r, (4.24) 

hold for all x € [a,h]. Note that it is enough to assume that the regularity 
condition holds for the constraints which are active at the local minimizer 
y. Indeed, suppose that I < r constrains, say Q^, . . . for simplicity, are 
active at the local minimizer y, and 



Sg\y;h^) 6gHy;h^) 

<5g2(y;hi) 6gHy;h^) 

Sg'{y;h^) 5g'{y;h^) 



SgHy;^] 
6gHy;y] 

sg'iy-M) 



7^0 



for (independent) h-' G D, j = Since the inequality constraints 

g^~^^,...,g^ are inactive, the conditions (|4.24p are trivially satisfied by 
taking A^+i = • • • = Ar = 0. On the other hand, since the inequality 
constraints g^, . . . ,g^ are active and satisfy a regularity condition at y, the 
conclusion that there exist constants Aj € M, j = 1, . . . , r, such that ()4.23p 
holds follow from Theorem 14. 4i Moreover, (j4.24p is trivially satisfied for 
j = l,...,/. 



5. Conclusion 

The FC is a mathematical area of a currently strong research, with nu- 
merous applications in physics and engineering. The fractional operators 
are non-local, therefore they are suitable for constructing models possess- 
ing memory effect. This gives several possible applications of the FCV in 
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describing non-local properties of physical systems in mechanics or electro- 
dynamics. In this note we extend the notions of Caputo fractional deriva- 
tive to the fractional derivative ^ D"'^ . We emphasize that this derivative 
allows to describe a more general class of variational problems and, as a 
particular case, we get the results of [2]. 

Knowing the importance and relevance of multiobjective problems of 
the calculus of variations in physics and engineering, our further research 
will continue towards multiobjective FCV. This is a completely open re- 
search area and will be addressed elsewhere. 
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